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Abstract 

We study the relation between the Roma and Zaragoza proposals for chiral 
fermions on the lattice. The fermion action in the Roma approach is shown to be 
equivalent to one of the Zaragoza type. This result is used to perform a mean-field 
study of the phase diagram for chiral Yukawa models based on the Roma action. 
The phase diagram is compared with the one based on the Zaragoza model with 
the most local choice for the fermion interactions. 


1 Introduction 


The formulation of a Chiral Gauge Theory (CGT) on the lattice suffers 
from the well-known doubling problem [lj]. Several proposals to deal with this 
problem have been reviewed in Refs. |2]||. In some of them one tries to apply 
the recipe successfully used in vector-like gauge theories: a Wilson-like term 
is introduced in an attempt to give a large mass to the unwanted doublers. 
In the Smit-Swift model for example [fij, a scalar held is introduced to write 
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down a gauge invariant Wilson-Yukawa term. When the scalar held acquires 
a non-vanishing vacuum expectation value, a mass term is generated for all 
the fermions. In the perturbative regime in the continuum, however, a fermion 
with a mass created by a Yukawa interaction does not decouple when this mass 
becomes large ]5|.|6j.|7[] . The same phenomenon appears to occur on the lattice. 
The electroweak S, U and A p parameters, for example, receive contributions 
from the doublers in the Smit-Swift model [Q. In the mirror fermion model 
0, supplementary, interacting physical fermions (the mirror fermions) are 
introduced. They make the model vector-like, with chiral properties at low 
energy, so that a conventional gauge invariant Wilson term can be used; the 
doublers get a mass of the order of the cut-off and the parameters are tuned in 
such a way that the masses of the mirror fermions are sufficiently high to have 
remained unobserved so far. These additional heavy fermions still contribute 
to the above-mentioned electroweak parameters, though ||. 


It seems that all models which try to produce heavy masses through the 
Higgs mechanism, to get rid of the unwanted fermions, are plagued by these 
non-decoupling effects. In the Zaragoza proposal 00-0 , the heavy fermions 
are avoided. Instead, the doublers are massless and non-interacting (see Sect. ^|). 
In the Roma approach |L3|, referred to as Roma I in Ref. j[[4|| , auxiliary 


fermions are used but the decoupling of these additional particles is achieved 
through a symmetry property [151. Both the Roma and Zaragoza approaches 
are gauge non-invariant regularizations. The Roma group has shown JT3[, how¬ 


ever, that one can recover a gauge invariant theory by an (in general) non- 
perturbative tuning of counterterms. The Roma and Zaragoza methods both 
preserve the global chiral symmetry and it has been shown that they reproduce 
the correct values for the fermionic contributions to the electroweak parame¬ 
ters S, U and (the “universal part” of) A p at one loop [§,[L6j 


The variant of 


the Roma model discussed in Ref. [[17]], known as Roma II ]F|], in which no 


redundant variables are used, and the Reduced Staggered Fermion Model of 
Smit and collaborators JT8|, should also correctly reproduce these parameters. 
However, as these models break the global chiral symmetry one would have to 
take into account the gauge non-invariant counterterms to achieve this 
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Recently, a lot of effort has furthermore been dedicated to proposals related 
to domain-wall fermions and fermions in the continuum [|T9|. However, the 
computation of the electroweak parameters S, U and A p using some imple¬ 
mentation of these ideas has not yet been considered. 


In this paperwe shall establish a direct relation between the Roma I and 
Zaragoza actions. This relation is exploited to determine the phase diagram 
of chiral Yukawa models based on the Roma I approach using mean-held tech¬ 
niques. This phase diagram is compared with the one for the Zaragoza model 
with the most local fermion interaction. 


For a concise presentation of part of the work described here, see Ref. [20]. 
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The study of such phase diagrams is important in order to investigate where 
a continuum held theory may be sensibly defined. Clearly, one would like to 
recover (the gaugeless limit of) the Standard Model as a continuum limit of a 
lattice model. Furthermore, it is not excluded, that certain properties of quan¬ 
tum held theories not anticipated on the basis of perturbative methods, like 
non-trivial hxed points, will provide a hint how to “improve” or extend the 
Standard Model. The search for such hints is another incentive to study phase 
diagrams of chiral Yukawa models. In addition, the answer to questions such 
as: what is the maximal fermion mass which can be generated by a Yukawa 
coupling, what is the effect of a strong Yukawa coupling on the upper and lower 
bounds on the Higgs mass, etc., will depend on where in the phase diagram 
of the, say, SU(2) x U(l) chiral Yukawa model one approaches the continuum 
limit. Most of these questions have already been studied for the Smit-Swift 
model |2]J and the mirror fermion model . In particular, in the strong cou¬ 
pling symmetric phase of the former model, some perturbatively unexpected 
physics was found |^3|. For a detailed discussion of possible continuum limits 
in general see the review of Ref. |24| . Some of these questions are currently also 
under investigation for the model with the local Zaragoza form factor. 

Here, we will only briefly discuss some properties which can be analysed in 
the context of our mean-held computations, such as which particles survive in 
the continuum limit. We would like to recall, however, that in the Roma and 
Zaragoza formulations the Standard Model is defined in the scaling region on 
the broken side of the small-y ferromagnetic-paramagnetic phase transition. 


The remainder of this paper is laid out as follows. In Section 2, the main 
characteristics of the Zaragoza proposal are recalled. Section 3 reviews the aux¬ 
iliary fermion method used in the Roma I action and establishes the equiv¬ 
alence with a model of the Zaragoza class. In Section 4 we calculate and 
compare phase diagrams of chiral Yukawa models in the mean-held approx¬ 
imation. Some subtle aspects of the Zaragoza phase diagram, extending the 
earlier work 0. are also discussed here. A summary and discussion can be 
found in Section 5. Details of some of the calculations have been collected in 
Appendices A and B. 


2 The Zaragoza proposal 


The philosophy of the Zaragoza approach Hail (see also Ref. |28|) to 
chiral fermions on the lattice is to tolerate the presence of the species doublers, 
but to prevent their communication with the real world: they are allowed to 
propagate but are kept massless and non-interacting. 
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The kinetic term for the fermions is taken to be the naive one, 

Skin = 51 ^( x ) K x,y^(y), (1) 

x,y 

where K is the usual lattice Dirac operator, 

F-x,y — X 55 TV i^y,x+p. ~~ $x,y+jj) ■ (2) 

It takes care of the propagation of all 2 d fermion species. 

The interaction terms are to be constructed in such a way that the doublers 
do not contribute. For this purpose, we introduce the quantity ipW defined in 
momentum space by 

fp (1) (p) = F(p)ip(p). (3) 


The “form factor” F(p) is required to be 1 for 
equals any of the doubler momenta. The idea is 
to the doublers at tree level. For the function F, 
choice was taken in Ref. 


Tlfl , for example) 


p = 0, and to vanish when p 
to prevent ip^ from coupling 
we usually choose (a different 


hp ) = n/w. 



0 G (— 7T, 7T], 


(4) 


see Fig. |T]. With this form factor (J|), the relation between the fields -0d) and 
ip is the most local one in x-space: ip^ l \x) sits in the center of the lattice 
hypercube at x and is equal to the average of the ip fields living on the corners 
surrounding it. This is easy to implement in a numerical simulation. 

The form factor F picks out the fermion ip^ at p = 0. Similarly, one can 
select the fermions ip^ centered around any of the doubler momenta pi = 0, 
P 2 = ( 7 r, 0,...), etc., by defining 

^\p) = F^(p)*P(p), F^(p)=F(p + Pj ) (j = 1,.., 2 d ). (5) 

The particular form factor of Eq. ([|) satisfies = 1, so that the kinetic 

term ip Kip reduces to J2j ip^ Kip^\ The fermionic measure in the path integral 
does not factorize, though. 

The interaction terms are then formed in a straightforward way, using ip^ 
instead of ip. For example, a lefthanded chiral gauge interaction is written as 
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Form factor Zaragoza local 




Fig. 1. Form factors F(jp\,p 2 ) on the Brillouin zone in two dimensions (momenta 
in units of tt), for the most local choice of Zaragoza form factor (]i]) (top) and for 
the Roma form factor (^), with parameter r = 1 for the Wilson term (bottom). 
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( 6 ) 


s ut = ^ - 1)*? 1 ! 1 + W 

X,fl 

^l\x + - 1 )^l\x) , 

where, as usual, ^4!^ = Pl^ 1 ^ with P L = |(1 — 75 ). Note that this interac¬ 
tion involves the lattice link variable Ul^ with the unit matrix subtracted. A 
Yukawa interaction between - 0 d) anc [ a scalar field <f> will be of the form 

= yYl (^l\x)^(x)^\x) +^\x)^ + (x)il; ( i ) (x)] . (7) 

X ' ' 


More complicated interactions can be constructed as well. The key point is 
to use only as a n interacting fermion field. In this way, all the vertices pick 
up form factors in momentum space which suppress the contributions coming 
from the doublers. One is left with only one physical, interacting fermion, 
while the doublers remain present but massless and decoupled. 


This decoupling can be shown to hold beyond tree-level as well ||12i|28i . One 
approach is to use the Reisz power counting theorem on the lattice [|29[ (see also 
Ref. |30|]), which formulates (sufficient) conditions for lattice Green’s functions 
to tend to the corresponding continuum Green’s functions (for undoubled 
fermions!) in the continuum limit. Anticipating a change of variables to be 
discussed shortly ([TSp , we write the denominator of the propagator (|T^) as 
(temporarily restoring factors of the lattice spacing a) 


D(ap ) = s 2 / (a 2 F 2 ) + m 2 F 2 


( 8 ) 


where we have defined 

s 2 = 51 sin 2 aP^, m = y($), (9) 


and 

F 2 = F + F 


( 10 ) 


(although we shall usually consider form factors which are real number valued 
functions). For application of Reisz’ theorem, sufficient conditions are that 
this denominator is a C°° function, tending to p 2 + m 2 in the limit a —► 0, and 
satisfying 

D(ap) > A (^2 ^ sin 2 (ap^/ 2) + m 2 j , (11) 


6 




on the entire Brillouin zone, for some constant A. (In addition, there are a few 
other, rather weak conditions on the numerators of the propagators and on 
the vertices which are satisfied in the present case and, in fact, in many other 
cases.) In the massless case, this inequality is easily verified for the local form 
factor F (|j): 


V a 2 F 2 


1 v^ 4 sin 2 (qp^/2) cos 2 (ap M /2) 
a 2 ^ Y\u cos 2 (ap u /2) 


> Ei s i n2 (W2). (12) 

[l a 


If rn 7 ^ 0, we need part of the s 2 /(2cl 2 F 2 ) term in (||) to yield the Am 2 term 
in (ED near the edge of the Brillouin zone, where the factor F 2 multiplying 
the m 2 -term becomes small. Writing, for example, 

D(ap) = s 2 /(2a 2 F 2 ) + [s 2 /(2a 2 F 2 ) + m 2 F 2 }, (13) 


the term between square brackets is larger than 


2 d m 2 


in the inner region 


where p M < 7 t/2 , and larger than (1 /am) 2 m 2 in the rest of the Brillouin zone, 
and one readily finds a value of A for which the inequality (|TTD is satisfied. The 
condition requiring C°° behaviour of this denominator is not satisfied for this 
local form factor, but it can presumably be replaced by a somewhat weaker 
condition which would allow application of the theorem in this case. With a 
less local form factor, for instance the one following from the Roma action 


see Sect. |3 . 1|) , the C°° condition can also be satisfied. 


In addition to this perturbative decoupling theorem, there are Golterman- 

0 


Petcher shift symmetries 
they take the form 


for the doublers [12,281. In momentum space 


V’(P) —> V’(P) + £$(P~Pi), 

4>{p) —* ip(p) + e5(p~Pi), (14) 

where e and e are independent parameters, i = 2 , ,., 2 d , and the pi are the 
doubler momenta. On the ‘doubler fields’ defined in Eq. (IJ), and their 
conjugates, these symmetry transformations (i = 2 ,.. ,2 d ) act as 


—> ip {3 \p) + eSijSpp - pj ), 

^ b \p) —* ^ b \p) + eSijS(p - pj ), (15) 


taking care of the decoupling of the doubler with label i (for more details, see 
Ref. p 8 |). Clearly, the prospective physical fermion ^M 1 ) j s no t decoupled. 


Evidently, the use of instead of ip in the interaction terms (£|[7|) does 

this means 


not respect the gauge invariance. As in the Roma I model [|13 
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that gauge fixing and ghosts have to be introduced explicitly in the lattice 
definition of the mo del. [3] (We are not considering gauge fixing and ghost 
terms nor gauge-variant counterterms here because they are not relevant for 
our purposes.) A global chiral invariance is preserved, though, because x ) 

is a linear combination of -0 fields at different sites. This global symmetry 
restricts the number of counterterms. For example, no fermion mass counter¬ 
terms are needed. 


As mentioned in the introduction, the Zaragoza method has been used to 
study the contribution of a doublet of fermions to the electroweak parameters 
S, U and A p on the lattice. At one loop the continuum result was repro¬ 
duced P~Ff1 , confirming the decoupling of the doublers. One should note that 
the gauge-variant counterterms, needed to recover the symmetries of the con¬ 
tinuum, do not contribute to these parameters at the one-loop level. A lattice 
simulation aiming at a non-perturbative demonstration of the decoupling and 
a non-perturbative determination of the leading fermionic contribution to the 
universal part of Ap is currently in progress 


It is interesting to note that a change of variables 

V>(p) -* F-'ipMp) 


(16) 


transforms the original action (B0®, schematically written as 


5 = 0AT0 + 0 (1) (17 - A')0 (1) + p0 ( 1 ) <F0 (1) 

= 0AT0 + i>F + (U - K)F 0 + p0A + <f>A0 (17) 

into 

S' = 0 (F + )- 1 AW 1 0 + 0 (f/ - A00 + p 0 <F 0 , (18) 


where U is the usual gauge covariant lattice Dirac operator. This transfor¬ 
mation, with constant Jacobian, shifts the form factors from the vertices to 
the fermion propagator. In Feynman graphs, this can be visualized as follows. 
Every fernrion-fermion-Higgs vertex from ([17]) contains a factor F 2 = F + F 
(and similarly for the interaction with the gauge field). These form factors 
can either be thought of as part of the vertex, as in (p~7|) , or be assigned to 
the fermion propagators connected to the vertex, as in fllS]) . In fact, in the 
formulation of Eq. (|T 8 |) the propagator (in the broken phase) has the form 


Sip) 


-is} + mF 2 2 
s 2 + m 2 A 4 


( 19 ) 


5 Ref. JrJ describes a proposal for a practical implementation. 






where 


4 = H 7m sin IV 

M 


( 20 ) 


For p —> 0, F —> 1 and the propagator corresponds to a massive fermion. 
Dne to the overall factor F 2 , however, the poles at the doubler momenta are 
suppressed. 


It is interesting to note that the “Zaragoza philosophy” mentioned at the 
beginning of this Section acquires a different interpretation in this formulation. 
While in the original form ( |T7| ) species doublers were present, albeit massless 
and non-interacting, they are absent from the -0-propagator following from 
Eq. (fE). This is not in disagreement with the no-go theorem [J1] because the 
action ([18]), in which the original 0H) is now the elementary fermion field, 
contains a non-local kinetic term. 


3 The Roma I approach 


One of the ideas to deal with the problem of chiral lattice fermions has been 
to introduce auxiliary righthanded partners Xn f° r the lefthanded physical 
fields -0L. Then a Wilson term can be formed which is used to decouple the 
species doublers. Similarly, an auxiliary lefthanded partner \l is introduced 
for each of the physical righthanded fields -0#, with a corresponding Wilson 
term. 

Ideas of this kind have been discussed in refs. |9j. 
uses. 

The mirror fermion model j9j . in which the auxiliary fermions are interacting 
fields, transforming under gauge transformations, is essentially vectorlike in 
nature; the auxiliary mirror fermions are physical. The idea is, however, to 
obtain an effective chiral theory at low energies by choosing the (both diagonal 
and off-diagonal) elements of the fermion mass matrix, or the matrix of Yukawa 
couplings to the Higgs field, in such a way that part of the spectrum resides 
at high enough energies. 

If the auxiliary fields do not transform under gauge transformations, the 
Wilson term breaks the symmetry. In the work of Smit |I2|^3|], however, pursu¬ 
ing a gauge-invariant formulation, the Higgs fields evoked by acting with gauge 
transformations turn the Wilson terms into Wilson-Yukawa terms, thereby 
restoring the gauge invariance. 


13,32,33], with different 
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In the Roma model [13], the auxiliary fields transform only under global 


transformations, so that the Wilson term breaks gauge invariance. The central 
idea of the Roma approach was, however, that the model has to be defined 
in the presence of gauge-fixing and ghost terms, so one must keep such terms 
also at the non-perturbative level. At the same time, one has to include all the 
counterterms needed to restore the BRST symmetry in the continuum limit. 
In this formulation, the Wilson term can be naturally maintained in its gauge 
non-invariant form. 


We will discuss this gauge non-invariant form of the action with auxiliary 
fermions, and will refer to it as the Roma I action. In the absence of gauge 
holds, the mirror model action takes the form of the Roma action if all the 
relevant couplings to the mirror fields are set to zero (he., the bare pure and 
mixing mass parameters and the Yukawa coupling for the mirror fermions). 


The fermion part of the Roma I action |T3|] consists of a minimally coupled 
kinetic term and a Yukawa interaction, supplemented with a Wilson term. It 
can be written in a condensed way as 


Sr = ^lUl^l + XrKXr + ^rUr^’r + X.l K X.l 

+ yif> L ® t l’R + y^R® + ^L 

+ ^l w Xr + ^rWXl + Xl w ^R + XrW^l- (21) 

The physical fermion fields are denoted by ip, their auxiliary chiral partners 
by X- is a ‘mixed’ Wilson term (with parameter r = 1), 

^Wx = bd x + P) + x( x — P) - 2x(x)] , (22) 

" V II 


the rest of the notation is the same as in Sect. We have not written the 
kinetic terms for the gauge, ghost and scalar fields, the scalar self-coupling, 
nor the gauge-variant counter-terms, because they are of no concern here. 

This action is invariant under the global chiral transformations 



— 5 ► ^L^L 

4>L ~ 


i’R - 

Qr i’R 

VY - 



U L , 

Ur, 




Xr ~ 

—> Xl L 

Xr 



Xr — 

Xr 




XL - 

—> n R 

Xl 



Xl ~ 

> Xl 


Xl L 






u„M;. 




Hr 

Ur, S2+ 





( 23 ) 
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All the terms in (fH| ) containing auxiliary fields are left ungauged. This 
feature distinguishes it from the mirror fermion model [!9j]. As a consequence, 
the local chiral symmetry is broken. On the other hand, one gains an invariance 
under a constant shift of the x fields, 


Xr( x ) —* Xr(x) + VR 
Xl(x) -> Xl(x) + VL 


from which the decoupling of these fields follows [|15 


(24) 


3.1 Roma a, la Zaragoza 


Now we make the following observation. As the auxiliary fields appear at 
most quadratically in <S R (H), it is possible to integrate them out exactly. 
This gives rise to an additional contribution 

- ^lWK-'W^l - ^rWK-'WtJjr (25) 

to the kinetic term of the 0 field, while the rest of the action remains un¬ 
changed. In the condensed notation of Eq. (|2TD , the action becomes 


Si, = iP l U l tJj l - i7 L WK' l W4> L 

+ ^rUrt/jr - 4>rWK- x W4>r 


+ ?/0 L <h0 fl + ?/0 i? <h + 0 L 
= 0(F+)- 1 J AF^ 1 0 + 0 (C/-X)0 

(26) 

+ T/0(<FPr + <f> + P L )0, 

(27) 

if we take 


f* = /Z7 

V s z + w z 

(28) 

(R for Roma), see Fig. [3]. with s 2 = s 2 (p) 
E m (1 - cos p^j. 

= )T /t sin 2 Pn and w — w(p) — 


The action 
given by 


is exactly of the Zaragoza form (18), with a form factor 


Note that Fr(p) satisfies the crucial form factor requirements Fr(p = 0 ) = 1 
and Fr(p = Pi) = 0 for the doubler momenta Pi. A difference between the 
form factors (|4|) and ( |28|) (cf. Fig. |lj) is that Fr — 0 only at the doubler 
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momenta whereas the local form factor F (f|) vanishes on the entire boundary 
hypersurface of the Brillouin zone. For small momenta p both form factors 
go as 1 — p 2 /8 , and in dimension d — 1 they are identical. Note that, for 
d > 2, F R in Eq. (p8|) is not well-behaved at the doubler momenta p*. It 
has conic behaviour, F R ( Pi — k) ~ Vk 2 , causing non-differentiability at these 
locations. This complicates an analytic study of the behaviour of F R in x- 
space, needed to determine the analogue of 'ip^(x) for this case. Numerical 
calculations appear to indicate that F R (x = Na,y = 0) and F R (x = y = Na ) 
fall off like (—1 ) N /N 3 in two dimensions. 

It is interesting to have a look at the the fermion propagator given by the 
action (^). At tree level, in the broken phase, one reads off 


Sip) 


—i4 + rriF R 
s 2 + w 2 + w?F‘ r 
- i£ + mF% 2 
s 2 + mV}) R ’ 


(29) 

(30) 


where m = y(4>). The Wilson term w 2 (p) in the denominator of Eq. (^9|) will 
decouple the doublers by the usual mechanism of giving them a mass of the 
order of the cutoff, while the conservation of global chiral symmetry is reflected 
by the absence of a ru-term in the numerator. This is to be contrasted with 
the usual Wilson fermions where such a term, although formally of order a, 
gives a finite contribution due to its inclusion in divergent loops. 


Alternatively, one may demonstrate the decoupling of the doublers by invok¬ 
ing the Reisz theorem or the Golterman-Petcher shift symmetry, as discussed 
in Sect. ^]. The denominator of the fermion propagator is now a smooth func¬ 
tion ( i.e ., it verifies the C'°°-condition of Reisz’ theorem), and in the symmetric 
phase the required inequality follows from 


s 2 /a 2 F 2 R = ( s 2 + w 2 )/a 2 = 


^sin 2 ap M + QT[1 — cos ap v ]Y 


= - E sin 2 (ap At /2) cos 2 (ap M /2) + (J2 sm 2 (ap u /2)f 


> — ^sin 2 (ap M /2) 


(31) 


In the broken phase, near the edge of the Brillouin zone, part of the w 2 -term 


in the denominator of the propagator (|29|) has to provide the mass term, as 
in the model with the local Zaragoza form factor, cf. Eq. (|T3|). In this outer 
region, where at least one of the momentum components is large, \p v \ > 7t/2, 
we have w 2 > 1 and the inequality follows easily. 
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The applicability of the GP shift symmetry is independent of the choice 
of form factor F, so we simply refer to the earlier remarks (14—15) for the 
decoupling of the doublers in Roma ‘a la Zaragoza’. 


4 Phase diagrams of Chiral Yukawa Models 


The equivalence of the Roma I action to a Zaragoza type of action can be 
exploited to study the phase diagram of chiral Yukawa models based on the 
Roma I action. 


The phase diagram for an SU(2)^ x SU(2)# Chiral Yukawa model in the 
Zaragoza approach with the local form factor F (f||) was studied in Ref. 
using both mean-held techniques and lattice simulations for rif = 2 fermion 
doublets. The mean held calculation, used in combination with a small or 
large-?/ expansion, was shown to give a rather good approximation to the 
Monte Carlo results as far as the structure of the phase diagram is concerned. 
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Here we shall study the phase diagram of the corresponding chiral Yukawa 
model with the Roma I action using the mean-held techniques of Ref. ]27| . 
Some details of this approach, which combines more standard mean-held tech¬ 
niques with expansions in the Yukawa coupling y, are described in Appendix 
A. We shall also present some new results for the Zaragoza phase diagram with 
the most local form factor, extending the previous analysis [|j]]. The results 
of the mean-held calculations for both phase diagrams are plotted in Fig. |^. 


The action for rif degenerate fermion doublets is given by 


k 

S = — — ^2 Tr ($ + (a;)$(a; + jx) + <h + (a; +/i) < h(a;)) 

" X,LL 

n f r _ _ 

+ + P) - A( x + P>hM x )) 

i= 1 Lx,n 

+ ^y^f\ x )(®( x ) p R + ® + ( x ) p L)^t\x) , 

x 


(32) 


for arbitrary choice of the form factor F used in the definition of the interacting 
fermion held The modulus of the Higgs held is frozen. This corresponds 
to infinite bare Higgs self-coupling A. With triviality in mind, we do not expect 
this to be a severe restriction. The phase diagram is given in terms of the two 
parameters y and k. Because of the symmetry $ —> —<f>, y —»• —y, the y values 
can be restricted to the range y > 0. For y = 0 we have the 0(4) model with 
critical points ±k c . The same is true for y = oo, as can be seen after rescaling 
the fermion helds, ip —> ip/y/y. There are some subtle aspects to the y —>• oo 
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Fig. 2. Mean-field phase diagram for SU(2) x SU(2) chiral Yukawa models with 2 
fermion doublets, for the local Zaragoza form factor (dashed curves) and for the 
Roma form factor (solid curves). Phases and multicritical points are indicated with 
capitals for the Roma case and with lowercase letters for the local Zaragoza case. 
The mean field result for the large-y FM-FI transition in the local case, which is 
reliable only for y g) 10 7 and would be invisible for the present choice of scale along 
the horizontal axis, has been extended to smaller y for clarity. The exact (mean-field) 
location of the FI-AFM transition in the region of strong Yukawa coupling, for the 
local form factor (dash-dotted line), is unknown. 

limit in the case of the local form factor, though; we will come back to this 
later. 

In the region of weak Yukawa coupling y, one finds paramagnetic (PM), 
ferromagnetic (FM) and antiferromagnetic (AFM) phases separated by the 
second order transition lines (cf. Appendix A) 



_ ?/ 2 n f AW) 

y 9 j P f > 



= ~l~ y 


2 n f j(W) 

o ] pa ? 


(33) 


where 
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7T 


AW) _ / CL P -T 

pf J (2vr) 4 s 2 (p) ’ 

— 7T 

r(w) = j ({4 P f 2 (p) F k(p) 
pa J (2vr) 4 s 2 (p) 


( 34 ) 


with F w (p) = F(pi + 7 r, • ■ ■ ,p 4 + 7 r). These curves meet in a point A with 
coordinates 


Va = 



k A 


-/r> / 4 /i->, 


(35) 


where 


/f> = /?> ± C' 1 


For the local Zaragoza form factor A one finds 

Av 

'pf 


/, ( 7 ) = 1.6179 10 " 2 , I^ ] = 8.4124 1(T 5 , 


hence 

y A = 7.8823/^/717, k A = -0.25261, 


while for the Roma form factor F, 


R 


AW) 

'pf 


2.5703 1(T 2 , 


I$P = 7.3343 10 ~ 3 , 


and 


y A = 7.3783/777, k A = -0.44964. 


(36) 


(37) 


(38) 


(39) 


(40) 


One sees that the local form factor F suppresses the contribution from the 
regions near the doubler momenta stronger than F r does: /(f) is very small 
for the local case fl39| ) and the transition line between the PM and AFM phases 
is almost horizontal. 

As mentioned in Sect. |3j. in the absence of gauge fields the mirror fermion ac¬ 
tion takes the form of the Roma action if all the relevant couplings to the mir¬ 
rors are set to zero. In Ref. the mirror fermion model with SU(2) x SU(2) 
chiral symmetry and rif — 2 was simulated in essentially^] this limit. Our 

6 A minor difference is that the simulation of Ref. was carried out for a small 
non-zero bare mixing mass = 1 — 8 K = 0(1/IV 2 ). 
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mean-field results agree quite well with Fig. 1. of Ref. |22|, showing the small- 
y region of the phase diagram, as can be seen after the corresponding mapping 
of parameters (y —> G^/2K with K ~ 0.125). 


Beyond the point A, one may consider the possibility of a ferrimagnetic (FI) 
phase, characterized by nonzero expectation values for the fields as well as the 
staggered fields (in the staggered average, fields at the odd sites contribute 
with a minus sign). The mean held method can also be applied to this situa¬ 
tion |27|, and the equations for the transition lines can be solved numerically. 
The reliability of the //-expansion involved, in this intermediate-?/ region, is 
monitored in terms of the quantity \y I- I n addition, one must 

check whether the free energy of the FI solution is actually lower than that of 
the competing FM and AFM solutions. 


For the local form factor, a FI phase was found [[27| , separated from the FM 
and AFM phases by second order phase transition lines. In Fig. ^], these lines 
have been plotted up to the point where \y I = 1- Comparison 

with data from the lattice Monte Carlo simulation of the model showed that 
the presence and location of the FI phase are predicted quite well by the mean 


held results [27 


We have used the same numerical mean held methods to look for a FI phase 
in the model with the Roma form factor. In this case, however, the assumption 
of second order phase transitions between FM and FI phases, and between 
AFM and FI phases, led to the inconsistent result of a FM-FI transition line 
lying below the AFM-FI line. 


In order to discover which of these candidate lines is (are) not acceptable, 
we made contour plots of the free energy, see Fig. [| Curves of constant free 
energy T < 0 were determined for FM and AFM candidate solutions. The 
FM contour curves run roughly parallel to the PM-FM transition and its 
continuation, where T = 0, and similarly the AFM curves parallel to the PM- 
AFM transition and its continuation. Subsequently, we determined the line 
of intersection points, at which the FM and AFM solutions have equal free 
energy. This line was found to lie in between the mutually inconsistent second 
order lines, implying that neither of the second order lines is acceptable: Above 
the line of intersection points of the free energy curves, no AFM solution is 
possible because the FM solution in the same point has lower free energy. 
Thus the AFM-FI line, found to lie above this intersection line, cannot be 
a transition between an AFM and a FI phase. Similarly, the FM-FI line is 
unacceptable. This is not yet sufficient to exclude a FI phase beginning at A. 
A FI phase bounded by first order transitions would still be possible. However, 
further inspection of the free energies for FM, AFM and FI mean held solutions 
in this region, using similar methods, confirmed that a FI phase is absent (at 
least close to the point A), while the FM and AFM phases are joined by a first 


16 







Fig. 3. Contour plot of the free energy for FM and AFM solutions in the region 
just beyond the weak-coupling PM phase. The solid lines are the mean-field results 
for the second order transitions between the PM phase, where the free energy is 
zero, and the FM and AFM phases. The mutually inconsistent (inverted) results 
for second order FM-FI and AFM-FI transitions are drawn as dashed lines. The 
dotted curves connect points at which FM or AFM solutions have a given free 
energy. The predominantly vertical curves are those for the FM solutions, with free 
energies 0 = /o > fi > ... > /s, decreasing from left to right. The predominantly 
horizontal curves correspond to the AFM solutions, with the same values fo ,...,5 for 
the free energy, now decreasing from top to bottom. Points at which the FM and 
AFM solutions have equal free energy are indicated with the symbol +. 


order transition. This transition line is determined by the line of intersection 
points in Fig. |3]. It remains possible, though, that this first order line has an 
end point with a FI phase opening up beyond it [Q. 


At this point it is appropriate to make a remark about the coexistence 
of different phases at the point A. It was stated in Ref. |P| that, in the 


mean-field approximation, the FM-PM transition line should continue beyond 
the point A as an AFM-FI transition, with the same slope, and similarly 
the AFM-PM line would turn into a FM-FI transition line with the same 
slope at A. The existence of the point A would thus automatically imply 
the presence of a FI phase. This is not correct, as was shown in Ref. |34[ . 
In fact, the FI phase may be absent, as in the model with the Roma form 
factor studied here, and if it is present then all four slopes at the point A 
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may be different, as in the model with the local Zaragoza form factor. The 
validity of these arguments does not depend on the details of the mean-held 
approximation, such as the expansion in y used to deal with the fermion 
determinant. Quantitatively, however, different values for the slopes will be 
found for different approximations to the fermion determinant. In particular, 
the conclusion that a FI phase is absent in the model with the Roma form 
factor might have to be revised. Apart from this, one should of course keep 
in mind the limitations of the mean-held approximation. Fig. 1. of Ref. |22[ 


reporting the earlier-mentioned Monte Carlo study of the mirror model in 
the limit that it is essentially equivalent to the Roma action, displays a FI 
phase beyond the point A, but the authors remark that this region was not 
investigated in detail. 


Similarly, the phase structure in the large-?/ region can be studied by combin¬ 
ing mean-held with a l/y expansion am (see also Appendix A). One expects 
PM, FM and AFM phases, as in the 0(4) model to which the model reduces 
for y —> oo. Assuming second order transitions between such phases, one finds 
transition lines given by the curves 


kf(y) 

k pj(y) 

where 


I _ ilVr(S) 

4 y 2 2 pf ’ 

_I _ llVr(S) 

4 y 2 2 pa ’ 


AS) j d 4 p s 2 (p ) 

pf J (27r) 4 F\p) ’ 

— TV 

t(S) _ j d4 P s 2 (p) 

pa J F 2 (p)F 2 (pY 


These lines meet in a point B parametrized by 
y B = \/n f I { . S \ k B = -l[ S) /4I ( _ S) , 


with 

t(S) _ t(S) i t(S) 

2 ± 2 p{ -'pa 


(41) 


(42) 


(43) 


(44) 


The integrands of /pf p a (f42j) are just the inverse of the integrands of 7^4 
Is is easy to show, however, that for any form factor for which the inte- 

> 0, so that the points A and B are 


grals converge, both I_ 




> 0 and 
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guaranteed to exist. (Here we have an amusing example of positive functions 
/, g satisfying both / / > / g and / 1/f > J l/g.) 


For the Roma form factor (|28|) , one finds 


4? = 348.01, 4f> = 158.70, 


' P f 


pa 


(45) 


hence 

y B = 13.759^7, k B = -0.66917. (46) 


As before, we also looked for a FI phase in the intermediate-?/ region by solving 
the corresponding mean held equations numerically. The result turned out 
to be analogous to that for the small-?/ equations: An anomalous inversion 
between the assumed second order FM-FI and AFM-FI transition lines, with 
free energy considerations deciding in favour of FM and AFM phases separated 
by a first order transition. Again, the mean-held scenario does not exclude the 
possibility that a FI phase opens up some distance away from the point B, 
while the conclusion about the absence of a FI phase at B may again be 
modified by an improved treatment of the fermion determinant. 


For the local Zaragoza form factor the integrals 4>fpa ([42]) are divergent, 
implying the absence of a PM phase at large (Unite) y in this approximation. 
FM and FI phases were found, though [^?J, separated by an almost horizontal 
transition line. The l/?/-expansion involved in this case was estimated to be 
reliable for extremely large y values only, y 4 10'. 


The absence of a PM phase is reminiscent of the model of Ref. [36]. There, 


the fermions at a site x are coupled to the Higgs fields at the vertices of the 
hypercube at x. In the present model, the Higgs field at a site x is coupled 
to fermions at the vertices of the hypercube at x, so one might be tempted 
by the similarity to assign this model to the same class. In the model of Ref. 
0 , however, there is no equivalence to an 0(4) model in the y —> oo limit, 
whereas in the Zaragoza model with the local form factor such an equivalence 
is established by rescaling the fermion fields, > ip^/y/y, after which the 

fermions decouple. On the basis of this equivalence one can argue that the 
present model should be of the “funnel” type |37j, with FM, PM and AFM 
phases in the large-?/ region. 


The resolution to this paradox is somewhat subtle but quite interesting. De¬ 
tails of the calculations, in mean field, are given in Appendix B. It turns out 
that the phase diagram of the model with the local form factor can be consid¬ 
ered as a limiting case of a funnel-shaped diagram, in which the PM phase has 
been pushed towards y — oo. It has not disappeared, but is contracted onto 
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the line y = oo. The line element (y = oo, —k c < k < k c ) is a second order 
phase transition line between the FM phase at large but finite y and the PM 
phase at y = oo, and the multicritical point B coincides with the critical point 
— k c of the 0(4) model. In particular, the values of the mean fields are found to 
decay with critical exponents (calculated explicitly in Appendix B), when this 
line element is approached along lines of constant k. It would be interesting 
to study this region, in particular the vicinity of the point (y — oo, k — k c ), 
by numerical simulation on the lattice. 

An interesting question is whether the AFM phase of the 0(4) model at 
y = oo extends into the finite-y region. The same numerical approach men¬ 
tioned earlier is applicable in principle to investigate this. In practice there is 
a technical difficulty, however. The method makes use of the values of certain 
integrals which diverge in the limit y —> oo. At the extremely large y values 
required by the reliability condition these integrals are very hard to evaluate. 
We have therefore not been able to establish the precise location of the AFM- 
FI transition with these methods, but there are indications that the transition 
exists (the dash-dotted line in Fig. |2]), with the AFM phase extending into the 
finite-y region. 

In conclusion, we have found that the SU(2) x SU(2) chiral Yukawa model 
(|32|) with Zaragoza fermions generically has a phase diagram of the funnel 
type. An example is obtained by taking the form factor (PB| ) describing the 
Roma action. If the form factor is chosen to be the most local one (f|), we get 
the limiting case in which the large-y PM phase is pushed towards the line 
y = oo. Alternative local form factors, for which the integrals (]42|) are also 
divergent, show similar behaviour. 


In general, the dependence of the phase transition lines on the form fac¬ 
tor can be sketched as follows. If the form factor F suppresses the doublers 
strongly, i.e., if F is relatively small in the vicinity of the doubler momenta, 
then the integral 1^ (|34[) is very small. The small-y PM-AFM transition 


line will then be almost horizontal. The integral 1^ ([I|), which receives its 
dominant contribution from momenta around the origin, will be affected only 
weakly, but still the small-y PM-FM transition will be curved relatively little 
in this case. The PM-FM and PM-AFM transitions in the strong coupling 
region, on the other hand, are determined by the integrals /pfp a (f42|), whose 


integrands are the inverses of those of /pf' pa . These transition lines will thus 
be curved relatively strongly. As a general trend, modifying the form factor 
causes the transition lines at weak coupling to bend downwards and those at 
strong coupling to turn upwards, or vice versa. 
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5 Discussion 


We have established a relation between the Roma I and Zaragoza methods 
to put chiral fermions on the lattice. Both proposals preserve the global chi¬ 
ral symmetry and have been shown to reproduce, at one loop, the expected 
values of the leading fermionic contribution to the electroweak parameters S, 
U and A p. The decoupling of the fermion doublers and the protection against 
mass counter-terms due to the global chiral symmetry make these approaches 
well-suited for a non perturbative study of Chiral Yukawa Theories and in 
particular of the scalar-fermion sector of the Standard Model. One can, for 
example, perform a lattice Monte Carlo study of the decoupling of the dou¬ 
blers, and also a non-perturbative lattice computation of the most significant 
contributions to Ap is possible within the fermion-scalar sector of the Stan¬ 
dard Model |25l , [26| . 

We have demonstrated that the action in the Roma I formulation (where 
auxiliary fermions are present) can be rewritten ‘a la Zaragoza’ with the form 
factor given in Eq. (^). In x-space, this form factor is spread out widely, con¬ 
trary to the most local possible choice usually made for the Zaragoza method 
(Eqs. (j3|-[f])) which is easy to implement in numerical simulations. 

An interesting question would be whether other globally chirally invariant 
fermion actions can be rewritten in the Zaragoza formulation, with some par¬ 
ticular form factor. 


The reformulation of the Roma I action as a Zaragoza model was exploited 
to calculate the phase diagram of chiral Yukawa models based on the Roma 
I action in the mean-field approximation (combined with expansions in the 
Yukawa coupling y or 1/y). The results agree quite well with the Monte Carlo 
simulations in the small-?/ region of Ref. (22)]. In general, we argued that (mean- 
field) phase diagrams of chiral Yukawa models with Zaragoza fermions are of 
the “funnel” type, with PM, FM and AFM phases in both the weak and 
strong Yukawa coupling regions. For a class of local form factors, for which 
the fermion interactions are within the range of a few lattice sites, our mean- 
field calculation suggests that the PM phase in the large-?/ region is contracted 
onto the y = oo axis, part of which is at the same time a phase transition line. 
One should keep in mind, however, that higher order corrections in the 1/d 
expansion might change this conclusion. 


Of course, in the scaling region the physical quantities are expected to be the 
same for any form factor which kills the doublers. In particular, the physics 
at the FM side of the strong FM-PM phase transition in the Roma case 
should be the same as in the model with the local form factor. As remarked 
in Appendices A and B, at strong Yukawa coupling y the fermions do not 
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survive the continuum limit, a similar conclusion as was drawn in the Srnit- 
Swift model. In the small-y scaling region, perturbation theory and mean held 
indicate that the fermion content is as expected, with only the physical chiral 
fermion surviving in the continuum limit. 

Let us conclude with a brief speculation about a different region in the phase 
diagram where one might contemplate taking a continuum limit, namely the 
scaling region on the AFM side of the AFM-PM phase transition. Although 
the AFM and FI phases could be lattice artefacts (absent for certain non-cubic 
lattice geometries), there appears to be no rigourous argument precluding the 
existence of a (possibly non-trivial) continuum limit at an AFM-PM phase 
transition [Q. If such a limit makes sense, it may be interesting to investigate, 
for example, a weak-coupling scenario with fermions which are naturally light 
on the scale of 246 GeV: the fermions, with a mass proportional to v in leading 
order would remain (almost) massless in, say, a mean-field approach, while the 
gauge boson mass squared, proportional to k{(fi+(j) x+ p,), would presumably not 
differ appreciably from its value in the weak FM scaling region. 
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A Appendix 


Here we present some details of the mean held computation. The most 

As 


significant characteristics of the approach have been described in Ref. [27 


our starting point we shall take Eq. 11 of that reference, giving the mean held 
free energy for the SU(2) x SU(2) chiral Yukawa model of Eq. 


T = —-^rlog Z = |(A 2 + A 2 ) —2 dk(v 2 — v 2 ) + av + av 

— ^[u(a + a)+u(a — a)] — C 0 I. (A.l) 

Here Go = 2' [d ^ 2 ^nf/2, rif is the number of fermion doublets, d is the dimension 
of space-time, and u(a) = log(2/ 1 (o;)/a) = a 2 /8 — ad/384 + 0(a 6 ) for SU(2), 
with Ii the modified Bessel function of order 1. T is a function of the mean 
helds a, v and A, collectively denoted as hi, and of the staggered mean helds 
a, v and A, denoted as hj , and is parametrized by the hopping parameter k 


22 





for the Higgs field and the Yukawa coupling y. The quantities a and v, and a 
and v, are mean field values for auxiliary fields needed to describe the scalar 
field, and A and A are related to an auxiliary field corresponding to the fermion 


condensate (see Ref. El for details). 


The function / comes from the fermionic determinant. It is obtained after 
a truncation of terms of higher order than quadratic in y (or 1 jy) in the 
mean held average of the Yukawa terms. In fact, the real expansion parameter 
here is ynf(i]f^ rather than y (and y~ 1 nf(fj^rj^) rather than r/ -1 in the 
large-?/ region, where y and y are auxiliary fermion fields [^7]]), so that this 


approximation should be valid in and near the paramagnetic phases ]27 ] . For 
illustration we shall consider the region of weak Yukawa coupling, where / has 
the form 


/ = 


d d p [s 2 (p) + ( m2 ~ m 2 )F 2 (p)F 2 (p)] 2 + m 2 s 2 [p)[F 2 (p) 


(2tt) 


log 


Kp) 


HereF(p) is the form factor, F w (p) = F(pi+tt, ... ,p d +7i), s 2 (p ) = sin 2 (p M ), 
and m = (m + + m_)/2 , m = (m + — m_)/2 with 


m ± = y 


^(cr ± a) — (A ± A) 



(A.3) 


The dots on the function u(a) denote derivatives. We have normalized Eq. 
(| A. 2| ) such that the free energy of the PM phase is zero. Concerning the mass 
nif (in lattice units) of the physical fermions in the scaling region of the FM 
phase, we can, as usual, distinguish a weak and a strong FM region. In fact, it 
can be shown from Eqs. ( [A. 1| . [A.2| ) that m 2 = m 2 + (= y 2 z 2 in the notation of 
Ref. [[27| ) in the small -y FM phase. For large y the result is rnj = y 2 / z 2 [^7 


To hnd the phase of the system for given y and fc, the free energy is min¬ 
imized with respect to the mean fields (for a recent concise discussion see 
Ref. HU), for which a necessary condition is that the saddle-point equations 


OF _ OF 

dhi dhj 


(A.4) 


are satisfied. A second order phase transition from a symmetric to a broken 
phase occurs when a negative mode develops and the mean fields acquire non¬ 
zero values. In general, such a transition is given by the condition 

det F"(h,h) = 0. (A.5) 


This equation is hard to solve for the case that there are broken phases on 


.(A.2) 
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both sides of the transition under consideration (for example, a transition be¬ 
tween antiferromagnetic (AM) and ferrimagnetic (FI) phases). Then Eq. ( A.5 ) 
has to be satisfied for the non-trivial mean field values given by Eq. (|A.4| ). If 
one wants to solve these equations perturbatively by expanding T in powers 
of the mean fields, then it is crucial |34|] to include terms up to at least fourth 
order in the fields. Alternatively, a numerical approach may be adopted to 
End the exact solution, as in Ref. [^7| and in the present paper. It should 
be emphasized, however, that the assumption of second order phase transi¬ 
tion lines should be checked afterwards by comparing free energies of different 
mean field values. 


If either side of the transition is PM, however, then Eq. (|A.5|) is taken at 
zero mean fields and can be solved easily. One can for example expand T to 
quadratic order in the fields, 


T = h+Mh + h+Mh + O(ih.h) 1 ) , 


(A.6) 


where 
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(A.8) 


with 


AW) _ 
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d d p F A (p) 
(2n) d s 2 (p ) 


T {W) _ j d d p F 2 (p)F 2 (p) 
pa J ( 2n) d s 2 (v ) ' 1 ’ 
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From Eqs. (|A.4|,|A.5|) one Ends the curves 
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as second order transition lines for the PM-FM and PM-AM transitions in 
the small-?/ region, respectively. For d — 4 they coincide with the formulas 
( |33| ) in the main text, and for y = 0 one recovers the phase structure of the 
0(4) model. 


The region of strong Yukawa coupling y can be studied similarly, with the 
integral / of Eq. (|A.2|) replaced by its large-?/ analogue (see Ref. |27|). One 
Ends PM, FM and AM phases separated by the lines 


k$\y) 

k { P l\y) 


! 1 AS) 

1 - 


r 


2 F/2] 


n f^ ) 




with 


AS) 

V 


d d p s 2 (p ) 
(27r) d F 4 (p) 


AS) = f dd P g2 b) 

pa J (2n) d F^p)F^p) ’ 


(A.ll) 


(A-12) 


provided these integrals converge. For the Roma form factor fl28|) they do con¬ 
verge in d — 4 dimensions. For the local Zaragoza form factor (ff|), however, 
they are divergent for d — 4, suggesting that large-?/ PM and AM phases 
are absent in this model [27]. In this case there are some subtle aspects to 
the phase structure in the limit y —> oo, however, which will be discussed 
separately in Appendix B. 


B Appendix 


As we have seen, the mean held determination of the phase structure of the 
SU(2) x SU(2) chiral Yukawa model with the local form factor ([|) 

F(p) = n 008^/2), (B.l) 

m=i 


in the strong Yukawa coupling region, is complicated by the divergence of 
certain integrals (|A.121) . It was concluded in Ref. [^7j] that there are no PM 
and AM phases at large y, in contrast with the 0(4) model to which the model 
should reduce in the limit y —> oo, cf. the discussion in Sect. 4. Here we analyse 
this issue carefully and resolve the paradox. 
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First we investigate, in the large-?/ region along the line k — 0, whether 
or not there is a FM solution to the saddle-point equations with lower free 
energy than the PM candidate solution given by A = v = a = 0. As we will 
see, the answer depends crucially on the form factor F. For the local form 
factor ( |B.1| ), there turns out to be a FM solution which is favoured relative to 
the PM solution for all large but finite y, as suggested by the (minus) infinite 
coefficient of the 1/y 2 terms in Eq. (|A.11|). 


We start from Eq. ( |A.1|) for the mean field free energy T and take the 
staggered fields equal to zero. The saddle-point equations to be satisfied by 
the mean fields are 


a 

v 

A 


4 dkv , 


u + A + 


A 2 ii 

~2FF 



d d p (u + A)s 2 

(2n) d (ii + A) 2 s 2 + y 2 F 4: ’ 


(B.2) 


where A = —Xyjil(a), and Co = 2^ d / 2 ^rif/2. 
For k = 0, 


F 



—7r 


log 


A 2 s 2 + y 2 F i 
y 2 F A 


(B.3) 


and Eq. ( |B.2| ) simplifies to 


a = 0 , 

v = A = —|A, 

1 = C 0 \l 0 (v/y ), 

y 

where the function 


Io(x) 


d d p s 2 

(27r) d x 2 s 2 + F 4 


(B.4) 


(B.5) 


is seen to increase with decreasing x. 


26 


















For not too large y, these equations ( |B.4| ) have a non-trivial solution, and 
using the inequality 


Vg > 0 : — l — < ln(l + q ), 

q + 1 


(B.6) 


with 


A V 
y 2 F i ’ 


(B.7) 


one can verify that any such solution has lower free energy T than the com¬ 
peting PM solution. This holds for both the Roma form factor F Ti (^B|) and 
the local form factor (|B.1| ). 


Whether there is a FM solution for arbitrarily large y depends on the small- 
x behaviour of the function ( |B.5| ). If it does not blow up for x —> 0, i.e., if 
J 0 (0) is a convergent integral, then I 0 (v/y) is bounded from above by / 0 (0) 
and the last equation of (|B.d| ) admits no non-trivial solution for large enough 
y, implying a transition to a PM phase at a finite value of y. This happens for 
the Roma form factor Fr (|28|) in four dimensions. 


For the local form factor (|B.1|) , however, the fact that Iq(x) blows up for 
x —> 0 means that for all y < oo there is a FM solution v to this equation. In 
fact, we will demonstrate that 


Io(z) 




0 ), 


(B.8) 


in dimensions d > 2 (and ~ x 1 for d = 1). This implies that the FM solution 
to Eq. ( |B.4| ) behaves like 

/ 1 \ 1/3 

v = A ~ f —j (log?/)*' 2 ^^ 2) (y —>■ oo) (B.9) 


(we shall assume d > 2 from now on), indicating an approach to a PM phase 
with a second order transition at y = oo, with critical exponents given by 

(H)- 

This argument can be generalized for the approach of the y —> oo limit along 
any line of constant k satisfying — k c < k < k c , where k c = 1/d is the 0(4) 
critical point. Since we are close to a PM phase, we can expand Eq. o to 
lowest order in the mean fields, provided k c — k is not too small. Using again 
the small-x behaviour (|B.8| ) of Iq(x), we find that 
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(B.10) 

(B.ll) 


for y —> cx) along lines of constant /c < fc c . The width of the small-(1///) region 
where this approximation is valid shrinks to zero for k —> k Cl and for k = k c 
we find instead: 


v 
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r \j 


r^j 


oi ~ 


v 3 ~ 



1/? (log ! ,) (2 '' 7 »<‘ i - 2 » , 

(B.12) 

3/7 

! (logy) (6/7)(rf - 2) • 

(B.13) 


When the point (y = oo, k = k c ) is approached from below along curves of 
the form l/y ~ ((/c c — k)/k c ) /3 (/? > 0), the mean fields go like 


v ~ a ~ (fc c - k)^~ 2)/3 , A ~ (k c -k)^ +1)/3 (/3 > 7/2), (B.14) 
v ~ a ~ (k c -kf /7 , A ~ (k c -k) 3l3/7 (0 </3 < 7/2), (B.15) 

up to logarithmic corrections. In the limiting case f3 —» oo, corresponding to 
approaching the point from below on the y = oo axis, one finds v ~ ce ~ A ~ 0, 
as expected for the 0(4) model in the PM phase. 

We see that in this scaling region the fermion mass rrij ~ y 2 /v 2 blows up 
when the phase transition is approached, so that the fermions do not survive 
the continuum limit. 


For k > k c one finds in the limit y oo along lines of constant k, for k — k c 
not too small: 


k — k r \ 1 


a 


A 


k c 

' k - kX 11 '(i \ 1 ' 2 


k c a 

V = TI ’ 


k r 


j 


log y\ 


k r 


k — kr 


d—2 


(B.16) 

(B.17) 


which reproduces the mean field results for the FM phase of the 0(4) model 
for y — oo. 

Note finally that in all cases considered here the asymptotic behaviour of 
the product of mean fields o:A 2 or vA 2 is ~ (l/y)| \og(y/v)\ 1[d ~ 2 \ as can also 
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be deduced directly from Eq. (11). 


In conclusion, we see that, for the local form factor (||), the vertical line 
element y = oo, k 6 [— k c , k c ] is a second order phase transition line connecting 
the large -y FM phase in the region y < oo to the PM phase of the 0(4) model 
at y — oo. This result also implies that the divergence of the integral 1)0 
( A.12|) becomes meaningless, because the transition line k^J (y) between PM 
and AM phases has shrunk to a single point B = (y = oo, k = —k c ) on the 
y = oo axis. 


The promised demonstration concerning the x —> 0 behaviour ( |B.8|) of the 
integral J 0 (x) for the local form factor ( |B.1| ) is as follows. We shall consider 
the function 


7T 

h(x) = x 2 I 0 (x ) = I 


x 2 s 2 


d d p 

( 2n) d x 2 s 2 + F 4 


x 2 s 2 


d d p 
( 2n) d x 2 s 2 + 


(B-18) 


For d — 1 this integral can be computed exactly. One finds 

T ( \ = j^p / cos 2 (p/2) y 1 = 2x 

'' J 2ti \ (4x 2 ) sin 2 (p/2) J 1 + 2x 


(B.19) 


For the general case d > 2 we shall show that, for small enough x, I\(x ) 
satisfies 


Oix 1/2 ^log^ < Ii(x) < C 2 £ 1/2 ^log 


d—2 


(B.20) 


for positive constants C i j2 . 

A simple order of magnitude estimate already hints at an x 1 ^ 2 behaviour. 
The dominant contribution to the integral I\(x) for small x comes from the 
regions where F A (5 x 2 s 2 , such that the integrand is 0(1). The x-dependent 
size of these regions determines how large this contribution is. In this way 
one finds that the dominant contribution is of order x 1//2 and comes from the 
regions where N (1 < N < d — 1) of the momentum components p^ are close 
to 7T, with the other p M away from both 0 and i r. 


Rigorous bounds on h(x) are obtained as follows. For the lower bound one 
has s 2 (p) > sin 2 p!, so that 


Ji(x) 


> 


7r 



—7r 


d d p 
(27 T) d 



cos 2 (pi/2) 
(4x 2 ) sin 2 (p!/2) 


II sm 4 (p /i /2) 

li =2 


-1 


(B.21) 


29 















The integration over p\ can be done using the d — 1 result (|B.19|) , and subse¬ 
quently the p 2 -integration gives 


I\{x 


> 


dp 3 
2vr 



2x 


2x + sin 2 (p 3 /2)... sin 2 (p rf /2) 


(B.22) 


For the upper bound, we divide the integration region into subregions Rj de¬ 
fined such that p G Rj if sin 2 pj > sin 2 p k for all k ^ j\ Sl^d d p/ (27r) d = Ylj fn :j ■ 
On Rj, s 2 < dsin 2 pj, and using symmetry and positivity of the integrand one 
obtains 



cos 2 (pj/2) 

(4 dx 2 ) sin 2 (pj/2) 


II sm 4 (p /i /2) 

^3 


-1 


< d 


7T 



— 7T 


d d p 

(2vr) d 



cos 2 (pi/2) 
(Adx 2 ) sin 2 (pi/2) 


II sm 4 ^/ 2 ) 

[i =2 


-1 




(B.23) 


leading to an expression similar to (|B.22|) . 

Using furthermore that (pk/'x) 2 < sin 2 (p fc /2) < (p k /2) 2 , we get 



/9 \ 2 ( rf - 2 )' 



G 

x(-) 

< h(x) < dG 

xVd 


\7 Tj 

J 


where 

GH = Jdk 3 ...}dk i ^^- 1? 

0 0 V J “ 


(B.24) 


(B.25) 


We see that the limiting behaviour of I\ (x) is determined by that of G[x]. 

In order to obtain an upperbound on G[x\ we divide the integration region 
into 2 d subregions Si by writing 

1 V2x 1 

J dk^ = J dk,j, + J dk^ (p — 3,..., d) . (B.26) 

0 0 ^ 
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On those regions S t where one or more of the k M are smaller then yj2x , we 
write 


I St - 


2x 


2x + k\ ... 


I 1 < a/2x . 


(B.27) 


The remaining contribution is 


1 1 /- 2 - 

]dk 3 ... Jdk d j 2x+ * k , 

\/2x y/2x 
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1 1 



y/2x y/2x 


V2x 

k 3 ...k d 



(B.28) 


To obtain a lower bound on G[x\, we write 


G[x\ 


i 


i 



(2 a ;)/ 3 (2 a ;)/ 3 


2x 


2x + H ... k?j ’ 


(B.29) 


guided by what we have seen is a momentum region providing a large contri¬ 
bution to the integral. Taking f3 = l/(2d — 4), we get 

G[l] 2 / / ^k^T, = ^ (yr^' 06 ;) ■ (B - 30) 

(2a:)/ 3 (2a:)/ 3 3 d \ \ ) J 


Combining the bounds ( |B.28| . |BT30| ) on G[x\ with the estimate (|B.2d| ) for 
h{x) we hnally arrive at the desired result ( |B.20| ) with 


/ 1 \ d ~ 2 n \ d ~ 2 

c, = . c 2 = ^ (-) . (B.31) 

For d = 2 it is not difficult to improve the lower bound to C\ = \[2 while 
for d = 3 stricter estimates lead to the slightly stronger bounds C\ = \/2/7r, 
C 2 = 3 5 / 4 \/2/7r. 
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